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■ ABSTRACT 

o ; 
oo , 
o ■ 

• We present infrared spectroscopy from the Spitzer Space Telescope at one disk 

■ position and two positions at a height of 1 kpc from the disk in the edge-on spiral 
j_| ' NGC 891, with the primary goal of studying halo ionization. Our main result 

■ — ' is that the [Ne III] /[Nell] ratio, which provides a measure of the hardness of 

the ionizing spectrum free from the major problems plaguing optical line ratios, 
is enhanced in the extraplanar pointings relative to the disk pointing. Using 
a 2D Monte Carlo-based photo-ionization code which accounts for the effects 
of radiation field hardening, we find that this trend cannot be reproduced by 
any plausible photo-ionization model, and that a secondary source of ionization 
must therefore operate in gaseous halos. We also present the first spectroscopic 
detections of extraplanar PAH features in an external normal galaxy. If they 
are in an exponential layer, very rough emission scale-heights of 330 — 530 pc 
are implied for the various features. Extinction may be non-negligible in the 
midplane and reduce these scale-heights significantly. There is little significant 
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variation in the relative emission from the various features between disk and 
extraplanar environment. Only the 17.4 fim feature is significantly enhanced 
in the extraplanar gas compared to the other features, possibly indicating a 
preference for larger PAHs in the halo. 

Subject headings: galaxies: ISM — galaxies: spiral — galaxies: individual(NGC 
891 — methods: numerical 



Introduction 



One of the most important changes in our way of thinking about interstellar gas in 
spiral galaxies in the past twenty years or so derives from the discovery that the ISM is 
much thicker than had been previously thought. Layers of great vertical extent (referred to 
as "extraplanar gas" or "gaseous h alos") have been found i n just about every com ponent of 
the ISM, especially in X-rays (e . g.. iTiillmann et al.ll2006bl : [Strickland et al.l 120041. HI (e.g., 

Oosterloo. Fraternali. &: Sancisil 120071). radio continuum (e.g.. iDahlem. Lisenfeld. &: Rossa 

2006 ), diffuse ionized ga s (DIG ; e.g.. lRossa fc DettmailboOSb : RandlllOod ). and dust ( Howk fc Savage 



19991 : lAlton et al 

flows 



le.g. 



19981. boooh. This gas may originate i n star-formation driven disk-halo 



Norman fc Ikeuchil Il989l ). prir aordial infall (IChiappini, Matteucci. fc Gratton 



19971 ) or accretion from satellite galaxies (Ivan der Hulst &: Sancisil l2005l ). Evidence for a 



connection to star formation comes from the correlation of the brightness and extent of 



formation, both within and among galaxies ( 


Rand 


1996. 


1998a 




Rossa & Dettmai 


2003a 


Tiillmann et al. 


2006a 




Strickland et al. 


2004; 


Dahlem et al. 


20061: 


Howk & Savage 


1999) 



On the other h and, some Galactic High Velocity Clouds have metallicities as low as 0.1 Zq 
(IWakkerll2004l ). which is not expected for gas originating in the disk. 



1.1. Emission Lines from DIG 



In the Milk y Way, the DIG (or Reynolds Layer) accounts for most of the ionized ISM 
(lReynoldslll990l ). and represents a significant power requirement. Energetically, only photo- 
i onization by rn assive stars can comfortably explain the intensity of optical line emission 
(lReynoldslll990l ). which raises the question of how an ionizing source concentrated in a thi n 
disk can be responsible for such thick ionized layers (iMiller fc Coxlll993l : lDove fc Shulllll994j ). 



Much information regarding the source (s) of ioniz ation have come from optical emission 



line ratios in the Milky Way and external galaxies (e.g.. lRandlll997l . ll998bl ; lGreenawalt. Walterbos. fc Braun 
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1997; Collins fc Rand 2001 



Otte. Gallagher, fc Revnoldsl 12002 



Haffner. Reynolds. &: Tufte 



I999I : iHoopes fc Walterboslboosi : iMadsen. Reynolds, fc Haffneilbood ). Studies of the spatial 
behavior of several diagnostic line ratios have supported the idea that photo-ionization by 
massive stars is the primary ionization mechanism, but crucial aspects of the data require a 
more complex picture. The most characteristic result is that the ratios of [S II]AA6716, 6731 
and [N II]AA6548, 6583 to Ha generally increase with distance from the suspected ionizing 
source, in the Milky Way and in external galaxy disks and halos. This was initially in- 
terpr eted as a fall in the ionization parameter , U, with distance from the ionizing source 



[e.g., iDomgorgen &: Mathij 1 1994 iRandl Il998bl ) but more recently has been attributed to 
an increase in the gas temperature in the DIG relative to th e immediate star forming envi- 



ronment due to a source of extra non- ionizing heating (e.g., iHaffner et al.lll999l : lOtte et al. 
2OO2I . using further information from [O II]A3727/Ha; in edge-ons), while even more recently 
the importance of radiation field hardening d uring propagation (due to the longer mean 
free p ath for harder photons) has been stressed (IHoopes fc Walterbosll2003l : IWood fc Mathis 
20041 . hereafter WM). 



Measurements of other line ratios present further difficulties for a simple photo-ionization 
picture. In several external spiral galaxies, both edge-on and more face-on, [O III]A5007/H/? 
is found to increase, or at least vary little, with distance from the ionizing massive stars (e.g.. 



Rand 



2OO3I : 



1998b 



200G 



Galarza. Walterbos. fc Braun]ll999l : lOtte et al.ll2002l : iHoopes fc Walterbos 



Miller fc Veilleuxll2003l ). An ir icrease has also been found in the Reynolds Layer for gas 
near the tangent point near I = 27° (iMadsen &: Reynoldsll2005l ). This behavior is unexpected 



for dilute radiation fields and has led to the consideration of secondary sources of ionization, 
such as shocks, in many of the above-referenced stud ies. Extra non-ionizing heating in halos 
may also be a factor, although lCoUins fc Randl (120011 ) find that this is an unlikely explanation 
i n NGC 5775 and UGC 10288. Exceptions to this b ehavior of [0III]/H/3 do exist, however 
jTiillmann fc Dettmai][200ol : iMiller fc Veilleuxll2003h . [O I]A6300/Ha, which depends on the 
neutral fraction of H through a charge exchange reaction, wh e n det e cted in the halos of edge 



ons is generally enhanced relative to th e disk (IRandl Il998bl. 



It is very weak in the Reynolds Layer (IRevnolds et al 



in one external non-edge-on spiral, M33 (IVoges fc Walterbol 



2OO0I : iMiller fc Veilleuxl [20oi ) . 



19981), an d has only been detected 
2006h . HeIA5876 is difficult to 



detect, but He I/Ha traces the relative ionization fractions of helium and hydr ogen and thus 



i s a u seful measure of the hardness of the radiation field. In the Milky Way, iMadsen et al. 



(120061 ) have detected the line in a few diverse environments, finding a range of helium ion- 
ization fractions, but in general the diffuse gas is in a lower i onization state than in HII 
regions. The line is dete cted in M33 with large uncertainties ( Hoopes fc Walterbos 20031: 
Voges fc Walterbos! I2OO6I ). while in the halo of NGC 891 helium is about 70% ionized (IRand 
1997h . 
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Spectroscopy of the DIG in irregular galaxies has led to a similar picture, with massive 
stars being the main source of ioniz ation, but shocks possibly c ontri buting as well [based on 
observations of 15 such galaxies by lHunter fc Gallagherl (119971 ) and iMartinI (Il997l )]. 



In summary, then, although radiation from massive stars dominates the ionization of 
DIG, many other factors may be at play: non-ionizing heating, shocks or other secondary 
ionization sources, and radiation field hardening during propagation. In addition, photo- 
ionization structure will depend on the how the ionization parameter changes with distance 
from the massive stars, and on the mean stellar spectrum emerging from the star forming 
regions. Further complications arise from uncertainties in the halo gas abundances, and also 
of course from dust extinction, which is not insignificant even at dist ances 1-2 kpc from the 



midp lane of several edge-on spirals, including particularly NGC 891 (IHowk &: Savagdll999 
2000l . and references therein). 



T he Infrared Spectro graph (IRS; iHouck et al.l |2004| ) on board the Spitzer Space Tele- 
scope ( IWerner et al.ll2004j ) opens a unique window on the ionization of diffuse gas through 
its access to infrared gas-phase lines, in particular the 12.81 /im [Nell] and 15.55 /im [Nelll] 
lines, and in this paper we present mid- infrared spectra of the disk and halo of NGC 891. 
The ratio of these lines provides a diagnostic of ionization that is relatively insensitive to 
extinction, gas-phase abundances, and temperature (being low excitation lines in warm gas) 
- three of the biggest sources of confusion for the optical lines. The first two ionization 
potentials (IPs) of Ne are 21.6 and 41.0 eV. With fewer complications, this ratio should 
provide excellent constraints on whether the diffuse gas is purely photo-ionized or requires 
an additional source by comparing against models. The utility of these lines as ionization 
diagnostics in Galactic and extragalactic HII re gions and starb u rsts w a s demonstrated usiri g 
the Infrared Space Obs ervatory (ISO) by, e.g., iThornley et al.l (120001 ) . iGiveon et al.l (120021 ). 
and IVerma et al.l (120031 ) . 



The aforementioned mod e ling b y WM uses the 3D Monte-Carlo radiative transfer code 
of IWood. Mathis. &: Ercoland (120041 ) to address several of the issues associated with photo- 
ionization. This code is able to employ more realistic geometries for the ambient gas - 
such as clumpy or fractal ISM distributions that may have a great impact on the ionization 
structure - than can be don e by the ID codes previously applied to the DIG problem (e.g., 
Domgorgen fc MathisI 1 19941 ). It also incorporates radiation field hardening, non- ionizing 
heating, and a range of possible input spectra and luminosities. WM find that spectral 



hardening reduces the need for the non-ionizing heating source invoked by, e.g.. lHaffner et al. 



(119991 ) ■ by producing higher gas temperatures and ratios of forbidden lines to Balmer lines, 
while suppressing He I emission. However, it is still not possible to match all of the line ratio 
data for the Reynolds layer and the halo of NGC 891 without invoking additional heating or 
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a secondary source of ionization. In this paper, we apply models tailored to the diffuse ISM 
density distribution of NGC 891 to constrain the sources of ionization using existing optical 
line ratios and the unique information provided by our Spitzer data. 



1.2. Emission from PAH Features 



Although our primary focus is on gas-phase lines, Spitzer spectroscopy also provides ac- 
cess to the mid-infrar ed spectral features that are generally attributed to pqlycyc l ic aromatic 



hydrocarbons (PAHs: iLeger fc Pueet 



1984 



AUamandola. Tielens. fc Barker 1985: Bakes. Tielens. fc Bauscl 



— — I n — - — ^ i — a — i — — i a — — — ■? ' = — 

2001al : iBakes et al.ll2001bl : iDraine fc Lill2007l ) . The identification of individual large molecules. 



or families of such molecules, responsible for the various features, as well as their creation, 
dest ruction, heating and emiss i on in different environments are all topics of current inter- 



est. iPeeters. Spoon, fc Tielend (12004 ) argue that PAH emission in general should be a good 
tracer of B stars. While PAH emission should be broadly connected to star formation, it is 
suppressed relati ve to the emission from larger grains in the immediate vicinity of recently 
formed stars (e.g, iLebouteiller et al.ll2007l ) . This is usually attributed to destruction of PAH 



molecules by harsh UV radiation fields. IWu et al.l (120061 ) find that the equivalent widths 
(EW) of several PAH features are lower in BCDs than in typical starburst galaxies, and 
find a weak anticorrelation of PAH EWs with [Ne III]/[Ne II]. A stronger anticorrelation is 
found with a parameter given by this ratio multiplied by the UV luminosity and divided 
by the metallicity. Thus, hardness, radiation field intensity and metallicity may all bear 
on PAH chemistry and emission. In galaxy halos the hardness, at least in the 20-40 eV 
range, is traced by the [Ne III]/[Ne II] ratio. The starlight intensity will clearly be weaker 
than in disks. As PAH emission is dominated by single-photon heating, their spectrum is 
not expected to va ry with starlight intensity for the kinds of environments relevant here 
dPraine fc LilboOTh . 



Further modification of the dust grain po pulation may oc cur through sputtering, grain 
collisi ons, and other processes (jPraind l2003l ). For instance, iJones. Tielens. fc HoUenbach 
(119961 ) find that shocks, with speeds as low as 50 km s~^, could be a source of PAHs through 
shattering of larger grains. Radiative acceleration of grains , with a size-depend ent efficiency, 
may also cause a modification of the grain siz e population ([Ferrara et al.lll991l) . with smaller 
grains staying closer to their formation sites (jPavies et al.lll998l ). Th ese effects could l ead to 
differences in the grain population between disks and halos [see also iHowk fc Savagd (119971 ) 
for a general discussion of processes that may lead to halo dust]. The vertical extent of PAH 
emission in halos and possible changes in the relative strengths of PAH features with height, 
and how these vary in galaxies with different star formation rates, could shed light on the 
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chemistry of these large molecules. 



Here we report spectroscopic detections of PAHs at 11.2, 12.0, 12.7, 16.5, and 17.4 /xm 
and a broad plateau from 16.5 to 17.5 fim in the halo of an external normal spiral, for the first 



time (the first spectroscopic PAH detections in any external halo were by lEngelbracht et al. 



(120061 ) for the starburst galaxy M82). PAHs emitting in this wavelength r ange are thought t o 
be generally neutral and larger than PAHs emitting shortward of 10 fim (jPraine &: Lill2007l ). 
Features at 11-14 /xm are thought to arise from C-H out-of-plane bending modes. The strong 
11.2 /im feature is believed to arise from compact, condensed PAHs, while the 12.7 Atm feature 
is thought to arise from more irregularly shaped molecules (Ivan Dishoeckl 120041). The red 
tail o f the 11.2 nm feature is thought to arise mainly from cations (iHudgins &: AUamandola 
19991 ). Ivan Diedenhoven et al.l (|2004| ) find that the shape of the 11.2 /xm feature, in terms of 
the asymmetry of the blue and red sides, and the strength of the smaller feature at 11.0 /im, 
varies with the type of exciting source (i.e., HII regions, Herbig AeBe stars, young stellar 
objects and planetary nebulae). 



Features in the 15-20 /im region have been discovered only recently. iMoutou et al. 



(120001 ) first reported the 16.5 /im feature, in Galactic sources, while t he 17.4 /im feature can 



be see n in the spectra of some of the Galactic sources observed by IVan Kerckhoven et al. 



(2000) and was also tentatively seen in the nuclei of M82 and NGC 253 bv ISturni et al. 
(|2000l). A broad p lateau from 15 to 20 /im was found by IVan Kerckhoven et al.l (120001 ). 



Peeters et al.l (120041 ) explore how these features change with the type of exciting source (i.e., 
HII regions, young stellar objects and evolved stars), and conclude that the discrete features 
in this wavelength range arise generally from prominent bands in neutral PAHs that should 
be larger than those responsible for the 10-15 fim features, while plateaus are due to blending 
of less prominent bands. They also fin d that the 17 . 4 p,m feature may arise from several 
candidate PAHs. In the dust model of iDraine fc Lil (l2007l ). the PAHs responsible for this 
feature are also very large (300-10^ C atoms). Er nission in this wavelength range from C- C-C 
bending modes of larger PAHs was predicted by I AUamandola. Tielens. &: Barker! (119891 ). 



As far as p r evious observations of PAH emission from galaxy halos are concerned, 
Irwin &: Madden] (120061 ) . using archival ISO data for the edge-on NGC 5907, detect mid- 
infrared narrow-band emission well into the halo. They argue that PAH features should 
dominate the emission in the filters used, and thus that the PAHs in this halo have a typical 
scale-height of 3.5 kpc, with emission extending up to 6.5 kpc. NGC 5907 is a galaxy with 
a low star formation rate and little evidence for a DIG halo (iRandl Il996l ). This result is 
therefore rather unexpected and needs to be confirmed. 

Finally, PAH features from 6.2 to 1 1 .2 nm in the disk of NGC 891 have already been 
detected by lMattila. Lehtinen. Sz Lemkd (119991 ). who find variation in the band ratios along 
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the major axis, suggesting changes in the dehydrogenation degree and charge of the PAHs. 

2. Observations 

The data were taken on 2006 February 2 and 2006 March 4 (program ID 20380; PI: 
R. Rand) using the staring mode of the IRS Short-High (SH) module on board Spitzer. 
A log of the observations is given in Table 1. The SH module is a cross-dispersed echelle 
spectrograph providing spectral coverage from 9.9 to 19.6 /im with a resolving power of 
~ 600. The aperture has dimensions of 4.7"xll.3". In staring mode the SH module nods 
between two pointings centered 1/3 and 2/3 of the way along the slit. 

Deep observations were made at locations in the halo of NGC 891 centered at a height 
of approximately z — ±1 kpc (hereafter "east halo" and "west halo" pointings) above a 
location in the disk at i? = 100" from the center of the galaxy. The reconstructed pointings 
indicate that the actual slit centers were at heights of 970 and 1120 pc on the east and west 
sides, respectively. A shorter observation of the disk location was also obtained. The exact 
coordinates of the pointings are given in Table 1, and the pointings are shown in Figure 1. 
The disk pointing falls between HII regions visible in Figure 1, although more embedded 
HII regions may contribute to the emission. The long integrations for the halo pointings 
meant that the data were divided into three Astronomical Observation Requests (AORs) 
per pointing, while the disk observation required only one AOR. Because of the bright 
infrared background, separate sky observations close in time to the target observations are 
necessary for IRS exposures of faint emission, particularly when the emission is expected 
to fill the slit. A deep sky exposure was taken on 2006 March 4 as part of this program. 
Unfortunately, the integration time is not as long as for the halo observations, resulting in 
the sky exposure dominating the noise in the halo spectra. No suitable sky exposure exists 
for the 2006 February 2 observations. However, only a small minority of the data were taken 
on this date (see Table 1) and since all lines of interest are detected with good sensitivity, 
these data were not used in the analysis. 
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Fig. 1. — Section of an Ha image of NGC 891 (IRand et al.lll99Cll ). Boxes show IRS pointings 
(for each pointing, two overlapping boxes, showing the two nods, are drawn). The three 
(yellow) boxes closest to the long-slit are the IRS SH pointings discussed here. The other 
(red) boxes in the disk show pointings for the SH and LH GTO data discussed. The sohd 
line shows the orientation of the slit for the optical emission line data discussed. The slit 
width is 2.25". The white blotch in the Ha image is an artifact of continuum subtraction of 
a foreground star. 
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The data were processed through version S13.2 of the IRS reduction pipehne, providing 
the Basic Cahbrated Data (BCD) products. Post-BCD processing was carried out according 
to the Infrared Spectrograph Data Handbook Version 2.0. Charge accumulation can be 
present in long integrations, evidenced by steady increases in the signal, but no significant 
effect was found in our data. The sky observation was subtracted from all target observations, 
and "rogue" pixels cleaned using the contributed IRSCLEAN_MASK software. For each of 
the AORs hsted in Table 1 and for each nod, the 2-D spectra were averaged using sigma- 
clipping at the 3(7 level. All remaining extreme pixels not removed by IRSCLEAN_MASK 
were replaced with averages of adjacent pixels. One-dimensional spectra were extracted 
using the package SPICE (version 1.3). Emission from NGC 891 is apparent over the entire 
aperture in every pointing, so the spectra were extracted and calibrated (with calibration 
tables, based on standard star observations, included in SPICE) using the full aperture 
width. Note that while version S13.2 of the pipeline allows extended source extraction, it 
is only fully accurate if the source uniformly fills the slit. The exact corrections for other 
emission distributions are impossible to estimate unless the source structure is known in 
advance. However, our comparison of our disk spectrum with nearby GTO pointings (see 
below) demonstrate that any such corrections will not affect the main results of this paper. 
A final spectrum for each pointing was formed by combining the resulting 1-D spectra for 
each AOR and nod. As the spectral orders overlap, the ends of each order were removed 
from the final spectra. The spectra are shown in Figure 2. 
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Fig. 2. — IRS SH spectra of the disk and east and west halo fields of NGC 891. Detected 
gas and dust phase features are indicated. 
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As we only had one disk pointing, we checked whether our gas-phase emission hne ratios 
were representative by examining several pointings near our disk position in the Spitzer 
archive from GTO data, taken on 2004 August 7 (program ID 97; PI: J. Houck). Both SH 
and LH (covering 18.7-37.2 fim at a resolving power of ~ 600 with an aperture of 11.1" 
x22.3") data were reduced. These pointings, spaced by about 11" in R.A. and 28" in Deck, 
are also shown in Figure 1. Integration times per nod are 31.5 and 14.7 s for the SH and LH 
modules, respectively. We refer to these pointings as A to E, with A being southernmost. No 
dedicated sky observations exist for these pointings, so we used sky observations from 2004 
August 8 (program ID 148). Because of the differing ecliptic latitude of these observations, 
the spectra had to be scaled to the expected sky level for the NGC 891 pointings using the 
predicted sky levels provided by the SPOT software before subtraction. Otherwise, the data 
were processed as above. 

Line intensities and EWs were measured with the IRAF0 program splot by summing 
pixels over the line extent and subtracting a linear baseline. Error bars on intensities are 
based on the noise in the spectra and the estimated uncertainties in the flux scale of IRS data, 
as described in the Infrared Spectrograph Data Handbook Version 2.0. Of the various sources 
of uncertainty in the flux scale, most affect the overall scaling independent of wavelength, 
while one is dependent on the position of the line in the order, and can be as much as 5%. 
We assume conservatively that it is everywhere at its maximum value. The total systematic 
error on intensities is then about 9%. The random noise in a single pixel is of order 10^^^ 
erg cm~^ s^^. For error bars on line ratios and EWs, we use only the random noise in the 
spectra and the wavelength-dependent error. 

We checked that the fluxes in our halo pointings cannot be due to disk emission convolved 
with the PSF of the IRS SH module using the Spitzer contributed software STINYTIM. The 
emission in the disk pointing is peaked around the midplane. At a representative wavelength 
of 17 /im, we estimate that a point source in the midplane would contribute 0.04% of its 
peak intensity in our halo apertures due to the wings of the PSF, whereas the halo emission 
we detect for all features is at least 3% of the disk emission. Therefore we expect very little 
contamination from the disk in the halo pointings. 



We also show optical emission line ratios in Figure 3 from the long-slit spectra of iRand 



(Il998bl ). The observations are described in that paper. 



^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Results 



3.1. Ionized Gas Phase Emission Lines 



In all three pointings, we detect [Nell] 12.81 /im, [Nelll] 15.56 /im, and [Sill] 18.71 
fim (Figure 2; Table 2). [Sill] has an IP of 23.3 eV, similar to [Nell]. We have also placed 
upper limits on [S IV] 10.51 fim emission. All of our disk line intensities fall between those 
for GTO pointings C and D (Table 3), as would naively be expected from the positions. 
To the extent that these limited measurements can be used to find the equivalent scale- 
height of an assumed exponential distribution of emission, the range of such scale-heights 
would be about 280 - 350 pc for [Nell] and [S 111], and 380 - 480 pc for [Nelll], given the 
error bars, and assuming no extinction (but see below). These change little if the average 
of the five intensities in the GTO data are used instead as a representative disk intensity. 
The [Nelll] /[Nell] ratio is 0.13 ± 0.01 in the disk and 0.31 ± 0.03 and 0.32 ± 0.03 in the 
east and west halo, respectively (see also Figure 3a). For the five GTO pointings A-E, the 
[Nelll] /[Nell] ratios are 0.12, 0.13, 0.14, 0.13 and <0.13, respectively. For [S lll]/[Ne 111], 
we find ratios of 2.1 ± 0.2 in the disk and 0.8 ± 0.1 and 0.9 ± 0.1 in the east and west halo 
respectively. For GTO pointings A-E, the ratios are 1.9, 1.9, 2.0, 2.8 and < 3.2, respectively. 
Hence, our disk results seem representative of the general trend in that part of the disk, while 
the two sides of the halo give similar ratios to each other. There is therefore a significant 
variation of these line ratios between the disk and the halo. Such variations we will refer 
to as halo-disk contrasts. There is no halo-disk contrast in [S 111] /[Nell], within the errors, 
consistent with their nearly equal IPs. The [Ne lll]/[Ne 11] results indicate a higher ionization 
state in the halo than in the disk. In §4 we explore whether this contrast can be explained 
by pure photo-ionization or whether a second source of ionization is required. We simply 
note here that the result cannot be explained by an additional non-ionizing source of heating 
in a pure photo-ionization model, as this ratio is temperature independent. 

In a perfectly edge-on, gas rich disk such as NGC 891, extinction in the midplane may 
be significant even at mid-lR wavel engths. To es t imate this effect, we cal culate total gas 



column densities from the CO map of lScoville et al.l (119931 ) and the HI map of lOosterloo et al. 



( 120071 ). We calculate H2 column densities from CO intensities ass u ming a conversion factor 



of 3 X 10^° mol cm ^ (K km s" 
Galactic relation between Ay and gas column density from 



as used by IScoville et al.l (119931 ) . We assume the 



Bohlin et al. 



Draine 



(Il978f ). convert Ai 



2003! ). and finally employ a 



to Ak using a standard Galactic reddening law (e.g., 
value of A]^i,m/AK of 0.4, co nsistent with determinations for many directions in the Galactic 
plane by iJiang et al.l (120061 ). Extinction in this part of the spectrum is very uncertain and 
depends on the strength of the 18 /im sili cate absorpt ion feature, which, while clear in the 
ISO spectrum of the Galactic Center from iLutzl (119991 ). is not evident in our disk spectrum. 
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Nevertheless, with these inputs we calculate about 2 magnitudes of extinction at 15 fim at 
our disk position. Extinction at a z = 1 kpc is negligible. The possible effect of reddening on 
[NeIII]/[Ne II] is eve n less well co nstrained. If we use the extinction at 12 fim towards the 
Galactic Center from iLutzl (119991 ). which exceeds that at 15 fim, the disk neon ratio would 
be further lowered. 

Both disk and halo [ Ne III]/[Ne II] ratios are typical of nearby starburst galaxies studied 
by iThornley et al.l (120001 ) . and much lower than the HII regions 30 Doradus and W51 these 
authors compare against - a contrast they attribute to aging effects. However, these authors 
find that DIG is not likely to contribute much to their global ratios for starbursts. With 
increasing distance to the starburst, the 14"x27" ISO SWS aperture should encompass a 
greater fraction of diffuse gas. This led to the expectation of a decrease in the line ratio with 
starburst distance under the assumption of a lower ratio for DIG, and such a decrease is 
not observed by them. However, if the ratio generally increases with distance from ionizing 
stars, as found here, but with still relatively low values in the DIG, then such a trend might 
not be expected. 
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Fig. 3.— Dependence of (a) [Ne III]/[Ne II], (b) [N II]A6583/Ha, (c) [S II]A6716/Ha, (d) 
[QI] A6300/Ha (e) [O II I]A500 7 /H/3, a nd (f) He I A5876/Ha on z. The optical line ratios are 
from iRandl f|l997n and iRand The Ha profile, normalized to unit intensity, from 



Randl (Il998bl ). is shown in (b). Horizontal error bars in (a) and (f) reflect the extent over 
which the data have been averaged. In the other panels, the intensities were averaged over 
317 pc. 
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3.2. PAH Emission Features 

We detect five PAH features, at 11.2, 12.0, 12.7, 16.5 and 17.4 /im, and a broad plateau 
from roughly 16.5 to 17.5 /im (Figure 2; Table 2). The halo-disk emission contrasts of these 
five features as ordered above are 0.068 ± 0.006, 0.04 ± 0.02, 0.048 ± 0.011, 0.056 ± 0.003, 
and 0.14 ± 0.03. Most features therefore have similar contrasts within the errors, consistent 
with the expectation that the PAH spectrum should not vary much with starlight intensity 
as long as it is not very strong. Only the 17.4 fim feature stands out as having a relatively 
high halo-disk contrast. The contrast for this line still stands out if the upper limits for the 
GTO data are used for the disk intensities (although upper limits on some other features 
complicate the comparison). It is also the only feature with a significantly higher EW in the 
halo (0.12 and 0.079 vs. 0.038), suggesting a larger vertical extent for this feature relative 
to the small, warm grains responsible for the continuum. Averaging the east and west 
halo measurements, the equivalent scale-heights of an assumed exponential distribution of 
emission would be about 330 — 530 pc for these features. For two magnitudes of extinction 
in the disk, the range becomes 210 — 270 pc. The scale-heights change little if we use instead 
the average of the disk intensities in the GTO observations (Table 3). For the 11.2 and 12.7 
fim features (the only ones detected in all five GTO pointings), they becom e 460 and 380 



pc, re spectively, assuming no extinction. We also note for comparison that iMattila et al. 



(1l999l ) found an intensity of the 11.2 fim feature of 2000 ± 350 erg cm~^ s~^ arcsec"^ for 
a 24"x24" field centered at our disk position using ISO, somewhat larger than our value of 
1400 ± 100 erg cm~^ s~^ arcsec"^, although the aperture sizes are not matched (again, the 
IRS SH aperture area is 53.11 arcsec^). 

CO emission at this location has a narrower vertical distribution. It is found to have a 



FWHM of about 225 pc in the high-resolution observations of IScoville et al.l (119931 ) . which 



would correspond to an exponential scale-height of 160 pc. Modeling of the HI disk by 



Oosterloo et al.l (120071 ) reveals a thin disk with a scale-height of <0.3 kpc in the inner disk, 
rising to 0.5 kpc in the outer regions, and a thick disk with 30% of the HI mass, with a scale- 
height increasing with radius from 1.25 kpc to 2.5 kpc. The PAH features may therefore 
have a vertical distribution comparable to the main HI layer, but the extinction correction 
would make their distribution intermediate between the CO and HI layers. 



Irwin fc Madded (120061 ) modeled the vertical distribution of emission in NGC 5907 in 
the ISO 6.7W filter with one to three components. It is impossible to say whether there 
is one or indeed more vertical components to the PAH distribution from our data, but 
nevertheless one can state that if the PAH distribution in NGC 891 is represented by a 
single exponential distribution, then the range of FWHMs for the features of 450-730 pc, 
assuming no extinction, is less than that of even the narrowest component (830 pc) in NGC 



- 16 - 



5907 (in a model with three Gaussian layers). With our data we of course cannot rule out 
an extended tail to the PAH emission distribution with a scale-height of a few kpc as found 
for NGC 5907 by Irwin & Madden, but such an extended halo seems very unlikely in the 
more actively star- forming NGC 891. 

Because four of the five PAH features show no significant halo-disk EW contrast, there 
is obviously no correlation with [Ne III]/[Ne II]. Note that we are only spanning a sma ll 
range, and a lower regime, of this ratio compared to the BCD sample of IWu et al.l (120061 ). 
where a weak correlation was found. However, the 17.4 fim EW is higher in the halo, 
where [Ne III]/[Ne II] is higher - the opposite of the weak trend found by Wu et al. for the 
6.2 nm and 11.2 fim features. The lack of EW halo-disk contrasts also puts limits on the 
modification of the grain size distribution by processes such as sputtering, shattering and 
radiative acceleration as discussed in §1. 

Also as discussed in §1, the origin of the 17.4 / xm feature is uncert ain - it could arise from 
several relatively large, neutral, candidate PAHs (Ivan Dishoeckll2004j ). The higher halo-disk 
contrast for this li ne may therefore indicate a halo favoring larger PAHs. Shattering of larger 
grains by shocks (jJones et al.lll996l ) in the halo is a potential source of such molecules. 



The spectral resolution of the IRS SH module permits an analysis of the shapes of 
PAH emission feature profiles. The profiles of the 11.2 /im feature, after subtraction of a 
linear continuum interpolated across regions adjacent to the line, are shown in Figure 4 (the 
two sides of the halo have been av eraged). The shape of this fea ture is found to vary with 
environment in the Milky Way by Ivan Diedenhoven et al.l (120041 . see their Figure 3). Most 
common is the Class An. 2 profile (relatively steep rise of the blue edge and fall of the red 
tail), typical of HII regions, non- isolated (i.e., associated with reflection nebulosity) Herbig 
AeBe stars and young stellar objects), while the Class Bn 2 (shallower rise of the blue edge, 
slightly redward shifted peak, shallower red tail, and the 11.0 /xm secondary feature shifted 
redward and broadened) is shown by many planetary nebulae and post-AGB stars. A few 
sources feature a hybrid class A(B)ii,2 profile (similar to An. 2 but with a shallower redward 
falloff and a slight redward shift of the 11.0 fim feature). In Figure 4, there is no discernible 
difference between the disk and halo profiles, suggesting little variation in the neutral/ionized 
mix or the excitation modes of the family of PAHs likely responsible for this feature. The 
profiles most resemble the An. 2 type, with a relatively steep blue rise, shallow redward falloff, 
and a clear secondary feature at 11.05 /im. The location of the peak of the secondary feature 
is more similar to that in the A(B)ii 2 profile. However, the main point is that there is no 
evidence for significant evolution of the PAHs responsible for the feature between the disk 
and halo. 
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Fig. 4. — The 11.2 iim feature for the disk (sohd hne) and the average of the east and west 
halo pointings (dotted hne). A hnear continuum has been subtracted from the spectra. 
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The 15-18 yum emission features a broad plateau on which the 16.4 and 17.4 /xm features 
sit. This part of the spectrum, after subtraction of a hnear continuum using regions adjacent 
to the plateau, is shown in Figure 5 for the disk and the average of the two halo pointings. 
The location of the plateau, and the strengths of the two features relat i ve to it and to each 
other, is much more reminiscent of the YSOs studied by iPeeters et al.l (120041 ) than the HII 
region in their sample. There is also a hint in the disk spectrum of the blue end of the second 
broad plateau seen in the spectrum of the ridge of H2 near the Herbig Be star LKHa234 by 
Peeters et al. from 18 to 20 /im, here seen as a shallow rise from 18 fi to 18.7 fim. Peeters 
et al. can reproduce both the HII region- and YSO-type profiles with mixes of laboratory 
PAH spectra, but the mixtures are not unique. 
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Fig. 5. — The 15-20 //m region for the disk (thick hne) and the average of the east and west 
halo pointings (thin hne) . A hnear continuum has been subtracted from the spectra. 
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3.3. Ho Emission 



H2 mid-IR emissio n lines arise from ro tational excitation in warm gas. It is usually at- 
tributed to PDRs (e.g., Ivan Dislioedal2004j ) but can also arise in shocks (e.g., iNeufeld et al. 



20061 ). Such emission ha s been mapped in the disk of NG C 891 using ISO Short-Wavelength 
Spectrometer (SWS) bv lValentiin fc van der Werj Jl999h . We detect the H2 S{1) J = 3-1 
17.03 fim rotationally excited line in the disk and halo. Our disk i ntensity is about 1.7 times 
higher than found for this line at the same position in the galaxy by lValentijn fc van der Werl 
(1l999l ). but the ISO SWS aperture size is not well matched to the IRS SH slit, being about 
seven times larger. Following the above analysis for the PAH features, if the emitting molec- 
ular hydrogen is in an exponential layer, then the measured halo-disk contrast suggests a 
scale- height of 390 — 460 pc. This result changes little if the GTO values are used. Again, ex- 
tinction may lower the scale-height significantly: using our extinction estimate in §3.1 would 
bring the scale-heights down to 220 — 260 pc, still slightly thicker than the CO-emitting 
layer. 



4. Ionization Modeling 

Figure 3 includes not on ly the [ Ne III] / [Ne II] ratio but also the optical line ratios along 
the long-slit (Figure 1) from iRandl ( 1998bl ). We wish to understand if these line ratios can 
be reasonably reproduced by photo-ionization from a population of massive stars residing 
in the disk. We focus on models that predict at least semi-quantitatively the observed neon 
ratio first, given the aforementioned advantages of this ratio as an ionization diagnostic, and 
then examine how well such models reproduce the optical line ratios. 

We calc ulated photo - ioniza tion models using the 3D Monte Carlo photo-ionization code 



described in lWood et al.l (120041 ). This code calculates the ionization and temperature struc- 
ture for 3D distributions of gas and ionizing sources for elements H, He, C, N, O, Ne, and 
S. We do not consider ions with ionization potentials above 54 eV (the second ionization 
potential of He), therefore we do not consider He^^ or radiation transfer of photons with 
energies above 54 eV. This simplification is justified by the lack of detection of He^^ emis- 
sion in the Reynolds layer and only its rare detection in HII regions. The models do not 
include the effects of dust on either the transfer of ionizing photons or the resulting line 
ratios. Observationally, from the CO and HI maps, modeling of the gas layers, and assump- 
tions about extinction in §3.1, Ay ~ 3 mag at 500 pc above a position 100" north of the 
galaxy center. Below this height, then, observed optical ratios may be somewhat affected 
by reddening. For instance, at 500 pc, the most reddened ratio. He I/Ha, may be reddened 
by about 0.7 mag. Other ratios are for lines much more closely separated in wavelength. 



- 21 - 



For the bulk of the optical emission line data points in Figure 3, then, reddening should be 
insignificant and will not affect our conclusions. The simulations can include heating over 
and above that from photo-ionization, according to any given prescription (see WM), but 
we have not included additional heating in these simulations. The elemental abundances we 
adopt by number relative to H are: H/He = 0.1, C/H = 140 ppm, N/H = 75 ppm, 0/H = 
319 ppm, Ne/H = 117 ppm, and S/H = 18.6 P2m^_Witli_the exception of S, these are the 
gas phase abundances in the local ISM used by iMathid (120001 ). The S abundance we use was 
found by WM to better match Milky Way and NGC 891 data. Note that the calculated S 
line ratios are highly uncertain due to the unknown dielectronic recombina tion rates . In ou r 
simulations we use the dielectronic recombination rates for S suggested by lAli et al.l (Il99ll ). 



All models have a central ionizing source and ionize a simulation volume that ex- 
tends to ±1 kpc in X, and z, with 65 cubical grid cells on a side. The ionizing spec- 
tra axe_40W£j_45kKj_aiid^ atmospheres from the WM-Basic library computed 
by [Sternberg. Hoffmann, fc PauldrachI (120031 ). The standard ionizing luminosity is 10^^ H- 
ionizing photons per second, except in one model where we consider a luminosity four times 
higher. Also, for one simulation, we use a " leaky" model in which the input spectrum is that 
from a spherical 40kK simulation which has its physical size set to allow 15% of the ionizing 
photons to escape. This leaky model suppresses the He-ionizing photons with energies above 
24.6 eV and hardens the spectrum in the H-ionizing continuum — see Figure 11 of WM and 
the accompanying discussion there. 



The density structure within the simulation grid is given by 

n{z) = 0.15exp[-|2|/0.66] + 0.04exp[-|z|/2.3] 



(1) 



where the densities are in cm~^ and distances are in kpc. This is a representative str ucture 
only, which we base on the HI data of Oosterloo et al. (2007) and lSwaters et al.l (119971 ). Our 
main results will not depend strongly on this assumed distribution. 

We note that the models predict [O HI] /Ha while our data are for [O III]/H/9. However, 
extinction should not affect the observed [OIII]/H/? since the lines are well matched in 
wavelength. Hence we expect trends with \z\ for [OIII]/Ha in the models to be comparable 
to those in [0III]/H/9 in the data. 
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Fig. 6. — Total projected relative intensity and line ratio maps for Ha, [NII]/Ha, [SII]/Ha 
and [Ne III] /[Nell] for a model with a 40kK ionizing source of luminosity 10^^ H-ionizing 
photons per second. 
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Fig. 7. — Vertical cuts at x — kpc (solid), 0.3 kpc (dotted), 0.6 kpc (dashed) and 0.9 kpc 
(dot-dashed) showing the variation of line ratios with z for the model in Figure 6. 
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We illustrate trends in the models by first showing in Figure 6 2D maps of the Ha 
intensity, [NII]/Hq!, [SII]/Hq!, and [NeIII]/[NeII] for our standard model with temperature 
40kK. Vertical intensity cuts that show line ratios at x offsets of kpc, 0.3 kpc, 0.6 kpc, and 
0.9 kpc are shown in Figure 7. Figure 6 shows that the [N II]/Ha and [S II]/Ha ratios increase 
with distance from the ionizing source. This is due to the progress ive hardening of the 
radiation field with increasing distance from the ionizing source (WM; iHoopes fc Walterbos 



20031 ). [Ne III] /[Nell] decreases with distance from the ionizing source due to the small 



number of high energy photons available at large distances to ionize Ne"*" to Ne^+ (only for 
the cut at X = 0.9 kpc do we see a hint of a rise in [Ne III]/[Ne II] with \z\, but it is clear 
from Figure 6 that this due to the idealized cone-like geometry of the simulations and cannot 
be considered a plausible match to the data; for the same reason, spikes and dips appear in 
most other line ratios at \z\ < 500 pc for x = 0.9 kpc, in clear disagreement with the data). 
[OIII]/Ha falls with \z\, Hel/Ha is overpredicted and [OI]/Ha is badly underpredicted for 
all cuts. The 50 kK model does produce [OIII]/Ha rising with \z\, but [Ne III]/[Ne II] ranges 
from 6 to 20 and falls with \z\, while the situation for He I/Ha and [O I] /Ha is not improved. 
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Fig. 8. — Vertical cuts at a; = kpc (solid), 0.3 kpc (dotted), 0.6 kpc (dashed) and 0.9 
kpc (dot-dashed) showing the variation of line ratios with z for the 40 kK leaky model with 
luminosity 10^^ H-ionizing photons per second. 
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The leaky 40kK model (Figure 8) shows falling [NeIII]/[NeII] with \z\ for all cuts. At 
X = 0.9, the values are at least comparable to the data. [OIII]/Ha and He I/Ha still fall 
with 1^1, again disagreeing with the data. Boosting the ionizing luminosity badly overpredicts 
[NeIII]/[NeII] and this ratio falls with \z\ between and 1 kpc, while all optical ratios are 
badly matched. The leaky, high-luminosity 40kK model best matches the run of He I/Ha 
(see discussion in WM), but again shows a neon ratio falling with \z\ and fails badly on the 
other optical ratios. In all the models considered here, [S III]/[Ne II] (not shown) is badly 
overpredicted. 

In summary, while the models can semi-quantitatively reproduce the behavior of [N II] /Ha 
and [SII]/Ha, they generally overpredict [Ne III]/[Ne II] and do not predict that this ratio 
should rise with \z\. Only at large distances in x from the ionizing source do the values begin 
to approach those observed (yet they still fall with \z\ in the model). The models also do not 
predict a rising [OIII]/Ha with \z\ as observed. The modeled [OI]/Ha values are generally 
too low, and He I/Ha is not well matched. These results demonstrate that hardening of the 
radiation field is insufficient to explain the data. An extra source of non-ionizing heating 
will also not affect the neon ratio as it is insensitive to gas temperature. 

One potential explanation for enhanced [Ne III]/[Ne II] in the halo would be the creation 
of higher levels of ionization due to increased collisional ionization. Such an enhancement 
could arise behind shocks in the halo. Although such shocks could not dominate the overall 
energetics of the diffuse ionized layers because of the large power requirements, they could 
provide the extra ionization at large heights above the disk. However, there are several 
difficulties with comparing our results to shock models. First, published models of steady- 
state shocks do not currently including predictions for infrared emission line intensities. 
Second, since these shocks are taking place in a low density medium, the length scales of 
the postshock cooling and recombination zones can be shown to approach the scale-height 
of the DIG. In such a ca se, it is necessary to consider the possibility of "truncated shocks" 



(iDraine &: McKed Il993l ) in which the integrated quantities are affected by the endpoint 
chosen for the recombination zone. Finally, given the hard radiation field produced by the 
postshock gas, the majority of the [Ne II] and [Ne HI] that forms occurs in the postshock 
recombination zone. Thus the column densities should be sensitive to the presence and 
spectrum of the external radiation field. A more complete investigation of the sensitivity 
and role of shocks in affecting the ionization structure at high latitude will be included in 
subsequent investigations. 
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5. Conclusions 

We have presented Spitzer IRS spectra of the disk and halo {z = ±1 kpc) of the edge-on 
spiral galaxy NGC 891. Our most important result is that the [Ne III] /[Nell] ratio is found 
to be enhanced in the halo relative to the disk, and this is very difficult to reproduce in pure 
photo-ionization models. 

Even this single observation presents a significant problem for models of DIG ionization 
featuring massive stars as the sole ionizing source. In a future paper we will explore shocks 
as a second source of ionization, but we remind the reader that other processes, such as 
cooling hot gas or turbulent mixing layers, may well be relevant and deserve further study. 

We have detected several PAH features in the disk and the halo. If in exponential layers, 
their scale-heights are in the rough range 330—530 pc, probably between the CO and HI scale- 
heights. However, extinction in the disk may reduce these scale-heights significantly. Most 
have similar halo-disk emission contrasts, suggesting little variation in the PAH population, 
as expected for low starlight intensities. There is also very little halo-disk contrast in EW 
for most of the features, putting limits on modification between disk and halo of the relative 
numbers of PAHs and the very small grains responsible for the continuum. Only the 17.4 fim 
feature stands out as having a higher halo-disk contrast and a higher EW than the others, 
possibly indicating a preference for larger PAHs in halos, although the origin of this feature 
is uncertain. The information on PAHs in halos is limited so far, but such work on the 
vertical variations of the various emission features has great potential for understanding how 
interstellar processes can transport grains and modify their populations. 

Our understanding of both the ionized gas and dust phases in halos will benefit from 
additional spectroscopic data in a range of halo environments, and we are in the process 
of acquiring such data with the IRS (program ID 40284; PI: R. Rand). In a future paper 
we will extend our NGC 891 data set with spectra at z = ±2 kpc, examine the halo of 
the more actively star-forming NGC 5775, where the rise of [OIII]/Ha with \z\ is even 
more problematic for pure photoionization models, and also the halo of NGC 3044, where 
a secondary ionization sourc e is not evidently required, as [OIII]/Ha clearly falls with \z\ 
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Table 1. NGC 891 IRS Observations^ 



AOR 


Date (2006) 


Pointing'' 


R.A. (J2000.0) 


Decl. (J2000.0) 


Integ. time per nod (sec) 


14595584 


March 4 


Disk 


2'^ 22™ 36.8= 


42° 22' 27" 


2x480 


14595840 


March 4 


Sky 


2^ 23™ 6.4= 


42° 22' 27" 


1x480 


14596096 


March 4 


Halo East 


2^ 22™ 38.6= 


42° 22' 22" 


9x480 


14596352 


March 4 


Halo East 


2^ 22™ 38.6= 


42° 22' 22" 


9x480 


14596608 


Feb 2° 


Halo East 


2^ 22™ 38.6= 


42° 22' 22" 


7x480 


14596864 


March 4 


Halo West 


2"^ 22™ 35.0= 


42° 22' 38" 


9x480 


14597120 


March 4 


Halo West 


2^ 22™ 35.0= 


42° 22' 38" 


9x480 


14597376 


March 4 


Halo West 


2^ 22™ 35.0= 


42° 22' 38" 


7x480 



^Spitzer program ID 20380 

''Pointing center of the field of view (each nod pointing is symmetrically offset from this position) 
'-'These data are not used 
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Table 2. Infrared Line Intensities and PAH Equivalent Widths 





Disk 


Disk 


Halo East 


Halo East 


Halo West 


Halo West 


Line 


Intensity* 


EW 


Intensity 


EW 


Intensity 


EW 






(/im) 




(/im) 




(/im) 


[SIV] 10.51 fim 


<1.9'^ 




<1.0 




<1.1 




[Nell] 12.81 urn 


320 ± 20 




12 ± 1 




11 ±1 




[Nelll] 15.56 urn 


41 ±3 




3.8 ±0.4 




3.4 ±0.3 




[SIII] 18.71 nm 


87 ±6 




3.0 ±0.3 




3.0 ±0.3 




H2 S{1) J = 3 - 1 17.03 nm 


49 ±3 




4.0 ±0.3 




4.3 ±0.3 




PAH" 11.2 


1400 ± 100 


0.85 ± 0.04 


105 ±9 


0.99 ± 0.05 


89 ± 7 


0.77 ±0.04 


PAH 12.0 Atm 


85 ±7 


0.043 ± 0.002 


2.5 ±0.9 


0.02 ± 0.01 


4.7 ±0.8 


0.04 ±0.01 


PAH 12.7 Aim 


600 ± 50 


0.32 ± 0.02 


36 ±3 


0.36 ± 0.02 


23 ±2 


0.21 ±0.01 


PAH 16.5 nm 


68 ±5 


0.08 ±0.01 


4.0 ±0.5 


0.08 ±0.01 


3.6 ±0.4 


0.08 ± 0.01 


PAH 17.4 fim 


32 ±2 


0.038 ± 0.002 


5.5 ±0.6 


0.12 ±0.01 


3.6 ±0.4 


0.08 ± 0.01 



''Units arc (10^-'^'' erg cm~^ arcsec"'^) 
''Upper limits arc 3(7. 

"^Intensities for the 11.2 fim feature include the secondary peak at about 11.0 fj,m. 
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Table 3. Infrared Line Intensities and PAH Equivalent Widths for GTO Data^ 



ABODE 



R.A. (J2000.0) 2^ 22™ 34.7= 2^ 22™ 35.5" 2^^ 22™ 36.5" 2^ 22™ 37.6= 2^^ 22™ 38.4= 



Decl. (J2000.0) 


42° 20' 57" 


42° 21' 50" 


42° 22' 19" 


42° 22' 46" 


42° 23' 14" 


Gas phase intensities'" 












[SIV] 10.51 fim 


<25° 


<25 


<25 


<25 


<25 


[Nell] 12.81 urn 


170 ± 14 


170 ± 14 


230 ± 18 


510 ± 40 


90 ±7 


[NcIII] 15.56 fim 


21 ±4 


23 ±4 


32 ±4 


66 ± 5 


<12 


[Sill] 18.71 /xm 


40 ±4 


43 ±4 


64 it 5.3 


180 ± 30 


30 ±4 


[SIII] 33.48 nm 


93 ±6 


150 ± 10 


170 ± 12 


280 ± 19 


81 ±6 


[Sill] 34.72 fim 


160 ± 12 


220 ± 15 


260 ± 20 


330 ± 20 


100 ±7 


H2 5(1) J = 3 - 1 17.03 ^J.m 


41 ±3 


36 ±3 


40 ±3 


51 ±4 


24 ±2 


PAHs intensities'^ and EWs (/xm) 












PAH'* 11.2 


940 ± 80 


750 ± 60 


1000 ± 80 


1300 ± 100 


740 ± 60 




1.04 ±0.05 


0.52 ± 0.03 


0.65 ± 0.03 


0.58 ± 0.03 


0.95 ± 0.05 


PAH 12.0 fira 


< 42 


< 42 


< 42 


< 42 


< 42 




< 0.05 


< 0.03 


< 0.03 


< 0.03 


< 0.05 


PAH 12.7 urn 


380 ± 30 


450 ± 40 


570 ± 45 


700 ± 60 


210 ± 17 




0.29 ± 0.01 


0.31 ± 0.02 


0.39 ±0.02 


0.36 ± 0.02 


0.245 ± 0.01 


PAH 16.5 /im 


< 19 


42 ±8 


53 ±8 


42 ±8 


< 19 




< 0.05 


0.07 ± 0.01 


0.09 ±0.01 


0.04 ±0.01 


< 0.05 


PAH 17.4 fim 


< 17 


< 17 


< 17 


< 17 


< 17 




< 0.05 


< 0.03 


< 0.03 


< 0.02 


< 0.06 



"Spitzer program ID 97 

''Units are (10~'^ erg cm~^ s~' axcsec"^) 

° Upper limits are 3(7. 

''intensities for the 11.2 /im feature include the secondary peak at about 11.0 /im. 




Right Ascension (J2000.0) 



